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Abstract 

Interpretation of measurements of the muon charge ratio in the TeV range 
! depends on the spectra of protons and neutrons in the primary cosmic ra- 

diation and on the inclusive cross sections for production of ir^ 1 and K in 
the atmosphere. Recent measurements of the spectra of cosmic-ray nuclei 
are used here to estimate separately the energy spectra of protons and neu- 
trons and hence to calculate the charge separated hadronic cascade in the 
atmosphere. From the corresponding production spectra of /i + and /i~ the 
ratio is calculated and compared to recent measurements. The com- 
parison leads to a determination of the relative contribution of kaons and 
pions. Implications for the spectra of and are discussed. 



1. Introduction 

The muon charge ratio in the TeV range has been measured by MI- 
NOS [H, 0] and more recently by OPERA [ij. Both analyses use an analytic 
approximation as a framework for making an inference about the separate 
contributions of the pion and kaon channels to the charge asymmetry. In 
this paper a more detailed derivation of the muon charge ratio is used for 
the analysis. The muon charge ratio is expressed in terms of the spectrum- 
weighted moments for production of n ± and K by protons and neutrons 
in the primary cosmic radiation, following the analysis of Lipari The 
analysis here accounts for the special contribution of associated production 
of charged, positive kaons. 
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This analysis also accounts for the effect of the energy dependence of the 
composition of the primary cosmic-ray nuclei. Measurements from ATIC j^] 
and CREAM [f| 0] indicate that the spectra of helium and heavier nuclei 
become somewhat harder than the spectrum of protons above several hundred 
GeV. This feature for helium was recently confirmed by PAMELA [§]. 

Because muon neutrinos are produced together with muons in the pro- 
cesses 

vr ± -> p^ + v^ipy) and ->■ p ± + z^(^), (1) 

these results also apply to and v^. In the TeV range and above the 
contribution of muon decay to the intensity of muon neutrinos is negligible. 
For reasons of kinematics, kaons are relatively more important for neutrinos 
at high energy than for muons. An additional goal of this paper is to draw 
attention to the implications of the muon results for atmospheric neutrinos 
in the TeV energy range and beyond. 



2. Muon charge ratio 

The excess of p + in atmospheric muons can be traced to the excess of 
protons over neutrons in the primary cosmic-ray beam coupled with the 
steepness of the cosmic-ray spectrum, which emphasizes the forward frag- 
mentation region in interactions of the incident cosmic-ray nucleons. The 
classic derivation of the muon charge ratio [9] considers muon production 
primarily through the channel p — >• 7r ± + anything. The atmospheric cascade 
equation for the intensity of nucleons as a function of slant depth X in the 
atmosphere is solved separately for N = n + p and A = p — n subject to 
the appropriate boundary conditions. For the total intensity of nucleons as 
a function of slant depth X (g/cm 2 ) 

ME) = MO) x exp(-^) (2) 

A N 

where the nucleon attenuation length is A N = A/(l — Z NN ) and A is the 
interaction length of nucleons in the atmosphere. The corresponding result 
for A(X) = p(X) - n(X) is 

A(X) = <Wiv(0) xexp(-^), (3) 



where 



p(0)-n(0) 1 1-Z pp + Z pn 1 

On = -^—^ 7 — r and - — = — - — . (4) 

p(0) + n(0) A_ 1 + Z pp + Z pn A N V ; 
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The Z-factors (like Z NN = Z pp + Z pn ) are spectrum- weighted moments 
of the inclusive cross sections for the corresponding hadronic process. For 
example, a particularly important moment for this paper is 

for the process 

p + air — > K + + A + anything. (6) 

The normalized inclusive cross section is weighted by x 1 where 7 is the in- 
tegral spectral index for a power- law spectrum and x = Ex/E p . Feynman 
scaling is assumed in these approximate formulas, so the parameters may 
vary slowly with energy, especially near threshold. However, the scaling 
approximation is relatively good because the moment weights the forward 
fragmentation region. 

2.1. Charged pion channel 

The next step is to solve the coupled equations for the production of 
charged pions by nucleons separately for I1 + (X) + U~(X) and for = 
U + (X) — U~(X). The solutions are then convolved with the probability per 
g/cm 2 for decay to obtain the corresponding production spectra of muons 
and neutrinos. The decay kinematic factors are 

-%1T ,'J: ~, (7) 



( 7 + l)(l-r 7r ) ' - cos0£ M ( 7 + 2)(l-rv) 
for muons and 

M -"- r " and £ - ^; r ^[ +1) (8) 



(7 + 1) ^ cos0E„ (7 + 2) 

for neutrinos. In each of Eqs. [7] and [8] the first expression is a low-energy 
limit and the second a high energy limit, where low and high are with respect 
to the critical energy e n . The ratio r n = m^/m 2 = 0.5731. The forms for 
two-body decay of charged kaons are the same with r K = 0.0458. 

The production spectra are then integrated over slant depth through the 
atmosphere to obtain the corresponding contributions to the lepton fluxes. 
Finally, the low and high-energy forms are combined into a single approxi- 
mate expression. 
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For example, for the flux of z/ M + the expression is 
<t> v {E v ) = MEu) 

I A 7TU 



B 7Tl/ cos(6)E 1/ /e 7T 1 + B Kv cos(6)E l/ /e K 

^■charm v 1 



1 + B chaiVm „cos(8)E l/ / e charm 



(9) 



Here 4>n{E v ) = dN/dln(E u ) is the primary spectrum of nucleons (N) evalu- 
ated at the energy of the neutrino. The three terms in brackets correspond 
to production from leptonic and semi-leptonic decays of pions, kaons and 
charmed hadrons respectively. The term for prompt neutrinos from decay of 
charm has been included in Eq. (see Ref. [10|) but will not be discussed 
further here. 

The numerator of each term of Eq. [9] has the form 

Z^i x BR lu x Zi v 

A iv — : ^ U<JJ 

with i = tt^, K, charm and BRi V is the branching ratio for i — > v. The first 
Z-factor in the numerator is the spectrum weighted moment of the cross 
section for a nucleon (N) to produce a secondary hadron i from a target 
nucleus in the atmosphere, defined as in Eq. [5j The second Z-factor is the 
corresponding moment of the decay distribution for i — > v + X, which is 
written explicitly in Eq. |8j The second term in each denominator is the ratio 



of the low-energy to the high-energy form of the decay distribution Uj. The 



forms for muons are the same, but the kinematic factors differ in a significant 
way (Eq. [7] instead of Eq. [8]). Explicitly, for neutrinos 



and for muons 



7 + 2\ /l-( r ^ +1 \ / A n -A 



B ^ U+iJ [i - VJy< +2 ) [khk/a-n)) ' (12) 

The forms for kaons are the same as functions of r% and A^. 

The separate solutions for n + — > n + + and 7r~ — >■ fj,~ + have the 
form 

rk (TP \± A (TP X 0-5(1 ±Q!,rffio) 

ME,) =ME,) 1+ mik ,^ (13) 
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where 



B 



± 

KfJ, 



B 



TTfX 



1 ± a n (35o 
1 ± c^a^do 



Here 



1 - Z 



i'/i 



J pn 



Z pp + Z 



0.909; & 



1 Z n + 7l + Z-x + Tt 



pn 



0.929; 



ou 



y 



0.165 



and 



1 - Ajy/A^ 



1 + 



Mflr/fi) 

ln(K/k N ) 



1.01. 



The numerical values are based on fixed target data in the energy range of 
hundreds of GeV [ill]. The factors differ by less than one per cent. To 
this accuracy, the charge ratio of muons can therefore be written in the form 



1 + /35 a m 



(14) 



where f n + = ( 
charged pions. 



1 - f35 a n 1 - f n + 
1 + /35r ) a 7r )/2 is the fraction of positive muons from decay of 



2.2. Leptons from decay of kaons 

The situation becomes more complex when the contribution from kaons 
is considered. In the first place, because the critical energies are significantly 
different for pions and kaons, the two contributions have to be followed sepa- 
rately. In addition the charge ratio of muons from decay of charged kaons is 
larger than that from pion decay because the process of associated produc- 
tion in Eq. [6] has no analog for forward production of K~ . Instead, associated 
production by neutrons leads to AK°. 

For the charge separated analysis of kaons it is useful to divide kaon 
production by nucleons into a part in which K + and K~ are produced equally 
by neutrons and by protons and another for associated production, which is 
treated separately. Then in the approximation that kaon production by pions 
in the cascade is neglected, the spectrum of negative muons from decay of 
K~ is 

; it? \- Z NK - A NK 



J NK 



+ B Kil cos(9)E^/e K 



(15) 
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There is an equal contribution of central production to positive kaons, 
but in addition there is the contribution from associated production. The 
total contribution of the kaon channel to positive muons is 



(En) A NK x 



Here 



and 



J pK+ 



1 + B+ cos(0)£ M / Cjc - 



J pK~ 



(16) 



B\ 



Z P K+ + ZpK- 

1 + (35 a K 



K » " 1 + AW(1 - ln((3)/ln(A K /A N )) ' 

Combining the expressions for fi + and yT from pions (Eq. IT3l) and from 
kaons (Eqs. IT5l and [TBI) . the muon charge ratio is 



cos(6)E ll /e 7T 

0--M 



+ 



+ 



\(l + a K f35 )A K JA^ 
1 + 5+ M C o^e)Eje K 

(Znk- I Znk) A Ktli j A n 



1 + B nfJl cos(6)E fl /e n 1 + B Kfl cos(6)Eje K _ 



(17) 



For the pion contribution, isospin symmetry allows the pion terms in the 
numerator and denominator to be expressed in terms of as defined after 
Eq. above. The kaon contribution does not have the same symmetry. 
Numerically, however, the differences are at the level of a few per cent, as 
discussed in the results section. 



3. Primary spectrum of nucleons 

What is relevant for calculating the inclusive spectrum of leptons in the 
atmosphere is the spectrum of nucleons per GeV/nucleon. This is because, to 
a good approximation, the production of pions and kaons occurs at the level 
of collisions between individual nucleons in the colliding nuclei. To obtain 
the composition from which the spectrum of nucleons can be derived we use 
the measurements of CREAM 0, 0] , grouping their measurements into the 
conventional five groups of nuclei, H, He, CNO, Mg-Si and Mn-Fe. 

Direct measurements of primary nuclei extend only to ~ 100 TeV total 
energy. Because we want to calculate spectra of muons and neutrinos up to 
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Figure 1: Left: three-population model of the cosmic-ray spectrum from Eq. 1211 compared 
to data 12h22|. The extra-galactic population in this model has a mixed composition. 
Right: Corresponding fluxes of nucleons compared to an E~ 2 ' 7 differential spectrum of 
nucleons and to the all nucleon flux implied by the Polygonato model (galactic component 
only) [H. 



a PeV, we need to extrapolate the direct measurements to high energy in 
a manner that is consistent with measurements of the all-particle spectrum 
by air shower experiments in the knee region (several PeV) and beyond, as 
illustrated in the left panel of Fig. [TJ To do this we adopt the proposal of 



Hillas [23] to assume three populations of cosmic rays. The first population 
can be associated with acceleration by supernova remnants, with the knee 
signaling the cutoff of this population. The second population is a higher- 
energy galactic component of unknown origin ("Component B"), while the 
highest energy population is assumed to be of extra-galactic origin. 



Following Peters [24] we assume throughout that the knee and other fea- 



tures of the primary spectrum depend on magnetic rigidity, 

R - g, (18) 

where Ze is the charge of a nucleus of total energy E tot = pc. The motivation 
is that both acceleration and propagation in models that involve collision- 
less diffusion in magnetized plasmas depend only on rigidity. The rigidity 
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Table 1: Cutoffs, integral spectral indices and normalizations constants dij for Eq. 1211 



determines the gyroradius of a particle in a given magnetic field B according 
to 

r L = R/B. (19) 

Peters pointed out that if there is a characteristic rigidity, R c above which 
a particular acceleration process reaches a limit (for example because the 
gyroradius is larger that the accelerator), then the feature will show up in 
total energy first for protons, then for helium and so forth for heavier nuclei 
according to 

E c tot = A x E NtC = Ze x R c . (20) 

Here En is energy per nucleon, A is atomic mass and Ze the nuclear charge. 
The first evidence for such a Peters cycle associated with the knee of the 
cosmic-ray spectrum comes from the unfolding analysis of measurements of 
the ratio of low-energy muons to electrons at the sea level with the KAS- 
CADE detector Q. 

In what follows we assume that each of the three components (J) contains 
all five groups of nuclei and cuts off exponentially at a characteristic rigidity 
R c j. Thus the all-particle spectrum is given by 

4>i(E) = E| =1 a it j E'^'i x exp 

The spectral indices for each group and the normalizations are given explicitly 
in Table [TJ The parameters for Population 1 are from Refs. p, 0] , which we 
assume can be extrapolated to a rigidity of 4 PV to describe the knee. In 
Eq. [21] 4>i is diV/ din E and 7$ is the integral spectral index. The subscript 
i = 1, 5 runs over the standard five groups (p, He, CNO, Mg-Si and Fe), and 
the all-particle spectrum is the sum of the five. 



E 



Zi Rc,j 



(21) 
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Figure 2: Solid line: charge ratio parameter So for the model with parameters of Table [T] 
Dashed line: same for Polygonato model [25| . 



The composite spectrum corresponding to Eq. [2TJ and Table [TJ is super- 
imposed on a collection of data in the left panel of Fig. [TJ No effects of 
propagation in the galaxy or through the microwave background have been 
included in this phenomenological model. For the two galactic components, 
however, a consistent interpretation could be obtained with source spectra 
7* ~ 1.3 for population 1 and 7* ~ 1.07 for population 2 together with an 
energy dependent diffusion coefficient D ~ E s with 5 = 0.33 for both com- 
ponents to give local spectra of 7 = 7* + 5 of ~ 1.63 and ~ 1.4 respectively. 
The extragalactic component comes in above the energy region of interest for 
this paper. We do not discuss it further here except to note that the last line 
of Table [1] gives the parameters for an extragalactic component of protons 
only. 

The spectrum of nucleons corresponding to Eq. [21] is given by 

<Pi, N {E N ) = A x 4>i(AE N ) (22) 

for each component and then summing over all five components. The nucleon 
spectrum is shown in the right panel of Fig. [TJ 

The energy- dependent charge ratio Sq(En) needed to calculate the muon 
charge ratio follows from Eq. [22] and Table [TJ To a good approximation, it is 
given by the fraction of free hydrogen in the spectrum of nucleons, as shown 
in Fig. [2] The fraction decreases slowly from its low energy value of 0.76 at 
10 GeV/nucleon [26| to a minimum of 0.63 at 300 TeV and then increases 
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somewhat at the knee. Note that, because of the relation among E tot , 
and R c in Eq. [2UJ the steepening at the knee occurs for nuclei at Z/A ~ ~ the 
energy per nucleon as compared to protons. Hence the free proton fraction 
rises again at the knee. 

Also shown for comparison in Fig. [2] by the broken line is the So param- 
eter for the rigidity-dependent version of the Polygonato model, which has 
a common change of slope A7 = 1.9 at the knee |25|. This gives rise to the 
sharp cutoff in the spectrum of nucleons for this model in the right panel 
of Fig. [TJ This version of the Polygonato model is meant to describe only 
the knee of the spectrum and the galactic component of the cosmic radia- 
tion. The behavior of the primary spectrum for En > 10 5 GeV/nucleon does 
not affect the charge ratio, which is measured only for < 10 4 GeV. It is 
therefore possible to consider the difference between the two versions of 5q in 
Fig. [2] as a systematic effect of the primary composition. 



4. Comparison with data 

We now wish to compare the calculation of Eq. [TTJto various sets of data 
using the energy- dependent primary spectrum of nucleons (Eq. 1221) with pa- 
rameters from Table [TJ There are two problems in doing so. First, expressions 
for the intensity of protons and neutrons from Eqs. [2] and [3] and the subse- 
quent equations are valid under the assumption of a power-law spectrum 
with an energy independent value of Sq. The assumption of a power law 
with integral spectral index of —1.7 is a reasonable approximation over the 
range of energies below the knee because it affects both charges in the same 
way. The proton-neutron difference, however, introduces an explicit energy- 
dependence into Eq. [171 that must be accounted for. We want to consider the 
energy range from 10 GeV to PeV over which the composition changes slowly 
with energy, as shown in Fig. [2J For estimates here we use the approximation 
5 (E N ) = <$ (10 x EJ. 

The other problem is that the data are obtained over a large range of 
zenith angles, and the charge ratio also depends on angle. The first MINOS 
publication |l| gives as a function of the energy of the muon at the 

surface. These data are shown in Fig . [3] along with older high energy data 



from the Park City Mine in Utah [27J and data at lower energy from L3 [28 



and CMS [29J. The figure shows three calculations of the muon charge ratio 



in the vertical direction that follow from Eq. [TJJ The highest curve assumes 



a constant composition fixed at its low energy value, 5q = 0.76 [26|. The 
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Figure 3: Muon charge ratio compared to data of CMS [29| and L3-C [28| below 1 TeV 
and to MINOS [l[ at higher energy. The L3-C data plotted here are averaged over 
0.9 < cos(0) < 1.0 for comparison with the calculation for vertical muons. See text 
for a description of the lines. 



middle curve is the result assuming the energy-dependent composition pa- 
rameter So(En) that corresponds to the parameterization of Tabled] (solid 
line in Fig. [2]), which is still higher than the data. Both the higher lines as- 
sume the nominal values of the spectrum weighted moments from Ref. fill ] . 
The lowest curve is obtained by reducing the level of associated production, 
by changing Z pK + from its nominal value of 0.0090 to 0.0079 

In order to look for the best fit it is necessary first to account for the 
dependence on zenith angle. The MINOS paper [l] does not give the mean 
zenith angle for each energy bin. However, because of the flat overburden at 
the Soudan mine where the MINOS far detector is located, there is a strong 
correlation between zenith angle and energy at the surface, as illustrated in 
Fig. 14 of Ref. [lj . Using this relation we estimate the effective zenith angle 
as a function of energy from cos(#) ~ 0.9 at 1 TeV to cos(6 l ) rs 0.5 at 7 TeV. 
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The points connected by lines in the multi-TeV range show the results of the 
calculation taking account of the dependence on zenith angle in Eq. [TT1 The 
lowest set of points has been adjusted to fit the MINOS data by varying the 
parameter Z pK +, which reflects p — > K + . 
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Figure 4: Muon charge ratio compared to data of CMS [H, OPERA 0] and MINOS 0. 
Measurements from the near detector of MINOS [3(J and from Park City [27j . 



The dependence of on zenith angle enters Eq. [T7] in the form 

E^cos(6). For this reason the muon charge ratio is often presented as a 
function of this combination. The data of OPERA are presented only in 
terms of the product E^cos(8). Because of the complex overburden at Gran 
Sasso, there is no simple relation between zenith angle and energy. The 
MINOS data are also presented in this form in Ref. [2J, but the mean energy 
for each value of cos(#) is not given. For a primary cosmic-ray composition 
that has no energy dependence, Eq. [T7| depends only on E^cos(9). The 
effect of the energy- dependence of composition can threrfore be assessed by 
comparing the calculation for various fixed values of 5q to the data, which is 
done in Fig. |H 
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The upper curve in Fig. [4] is the same as the corresponding curve in Fig. EJ 
plotted for a constant composition with 5 = 0.76, its value at 10 GeV. The 
parameter Sq decreases from 0.71 at 100 GeV/nucleon to 0.68 at a TeV, and 
from 0.64 at 10 TeV/nucleon to less than 0.62 at 100 TeV. The full curve 
through the data in Fig. [4] is evaluated for 5q = 0.665. The two broken lines 
in the low energy region are plotted for 0.71 and 0.69, while those at high 
energy are for 0.64 and 0.62. A more precise comparison between data and 
calculation could be made given complete information about the distribution 
of energy within each bin of zenith angle, but it is clear that the data from 
the various experiments are reasonably consistent with each other and with 
the present calculation. 



5. Summary 

The muon charge ratio is sensitive both to the proton excess in the spec- 
trum of primary cosmic-ray nucleons and to the value of Z p x+- Using recent 
data on primary composition, we find a proton excess that decreases steadily 
from 10 GeV/nucleon to 500 TeV. This portion of the cosmic-ray spectrum 
produces muons from a few GeV to well over 10 TeV. Assuming associated 
production (Eq. [6]) to be the major uncertainty, a level of associated produc- 
tion in the range Z pK + = 0.0079 ± .0002 is required to fit the observed charge 



ratio. For comparison, the nominal value [ll( is Z p x+ = 0.0090. Keeping the 



nominal values of all other parameters, the fit here corresponds to a ratio 

= Z vK , + Z vK - 0.0079 + 0.0028 
' Z pn+ + Z pT - 0.046 + 0.033 V ; 

It is interesting that analyses of seasonal variations of TeV muons by 



MINOS [31| and IceCube [32] also suggest a somewhat lower value of Rk-k 
than its nominal value of 0.149. On the other hand, the value in Eq. [23 
still represents a significant contribution from the K + decay channel. If the 
energy- dependent composition of the Polygonato model is used instead, a 
good fit is obtained with Z pK + = 0.0074, which reflects the somewhat higher 
value of 5q in the relevant energy range (0.68 as compared to 0.64). The 
fraction of kaons would be correspondingly lower (Rk/tt = 0.129). 
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In the analysis of MINOS, and also in that of OPERA, the muon charge 
ratio is written in the form 



1 + B Vfl cos(9)E fl /e 7T 1 + B Kfl cos(9)E^/e K _ 

, (1 - L+) | (1- f K+ )A Kll /A nlI 

1 + B nit cos(9)E lx /e n 1 + B Kfl cos(9)E^/e K 

with Axn/A nfl = 0.054. The more correct Eq. [T7|has a different form for the 
contribution from kaons. In the MINOS analysis the fitted values of the two 
free parameters are /+ = 0.55 and f£ = 0.67. For E^ ~ TeV, 5o ~ 0.64 for 
primary energy per nucleon of 10 TeV. Thus /+ = |(1 + a n (35o) = 0.55 in 
agreement with the MINOS analysis. A precise comparison with the MINOS 
value for f^- is not possible for the reasons explained after Eq. [T71 However, 
numerical differences are at the level of a few per cent. For example, p=j 
0.95 x B Kil . From Eqs. Hiand [161 the value of f K + = <M£ M ) + / '(<M^) + + 
(f)K{E^)~) ~ 0.69 in the TeV region for 5q ~ 0.64. However, if the expression 
for the kaon contribution in Eq. [T7] is expressed in terms of f£ there is an 
additional multiplicative factor less than unity. Thus, although the forms are 
different, the fits are much the same. 

The role of kaons is relatively more important for neutrinos than for 
muons. Because the muon mass is close to that of the pion, the muon car- 
ries most of the energy of the decaying pion. Kaons split the energy almost 
equally on average between the \i and the v^. The steep spectrum enhances 
the effect so that kaons are the dominant source of muon neutrinos above a 
few hundred GeV. Forward production of K + is therefore particularly im- 
portant. The effect is illustrated in Fig. [U in which the ratio v^jv^ is plotted 
for the same sequence of assumptions as in the plot of the muon charge ratio 
(Fig. [3]). The implications of the muon charge ratio for neutrinos will be the 
subject of a separate paper. 
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Figure 5: ratio of v^/v^ calculated with the same parameters as for the charge ration of 
muons. The dashed line shows the results with parameters of the Polygonato model 25 1. 
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